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Abstract 

A shape-function independent relation is derived between the partial B — * X u de- 
cay rate with a cut on P + = Ex — \Px\ < A and a weighted integral over the nor- 
malized B — » X s j photon-energy spectrum. The leading-power contribution to the 
weight function is calculated at next-to-next-to-leading order in renormalization-group 
improved perturbation theory, including exact two-loop matching corrections at the scale 
Mi ~ a/^Aqcd- The overall normalization of the weight function is obtained up to yet 
unknown corrections of order a1(m^). Power corrections from phase-space factors are 
included exactly, while the remaining subleading contributions are included at first order 
in Aqcd/ 77 ^- At this level unavoidable hadronic uncertainties enter, which are estimated 
in a conservative way. The combined theoretical accuracy in the extraction of \ V u b\ is at 
the level of 5% if a value of A near the charm threshold can be achieved experimentally. 



1 Introduction 



In the recent past, much progress has been made in the theoretical understanding of inclusive 
charmless B decays near the kinematic endpoint of small P + = E x — \Px\, where E x and Px 
are the energy and momentum of the final-state hadronic system in the £?-meson rest frame. 
In B — > Xg'-f decays the P + variable is related to the 5-meson mass and the photon energy, 
P + = Mb — 2E^, and the measurement of its spectrum leads directly to the extraction of the 
leading hadronic structure function, called the shape function [1, 2, 3]. The P + spectrum in 
semileptonic B — > X u l~u decays, on the other hand, enables us to determine \V u b\ [4, 5, 6], 
but this requires a precise knowledge of the shape function. One approach for measuring \ V u b\ 
is to first extract the shape function from the B — > X s 7 photon spectrum, and then to use this 
information for predictions of event distributions in B — > X u l~v. A comprehensive description 
of this program has been presented in [7]. Equivalently, it is possible to eliminate the shape 
function in B — > X u l~u decay rates in favor of the B — * X s ^ photon-energy spectrum. This 
idea was first put forward in [2] and later refined in [8, 9, 10, 11]. Partial B — > X u l~v decay 
rates are then given as weighted integrals over the B — > X s ^ photon-energy spectrum, 

exp. input theory exp. input 

where the weight function W(A, P + ) is perturbatively calculable at leading power in 
A comparison of both sides of the equation determines the CKM matrix element \V u i\ directly. 
For the measurement of the left-hand side to be free of charm background, A must be less 
than Mp/M B m 0.66 GeV. However, the P + spectrum in B — > X u l~v decays displays many of 
the features of the charged-lepton energy spectrum, so that it is not inconceivable that the cut 
can be further relaxed for the same reasons that experimenters are able to relax the lepton cut 
beyond the charm threshold. We stress that for an application of relation (1) a measurement 
of the B — > X s 7 photon spectrum is needed only for E 7 > \(M B — A) 2.3 GeV (or slightly 
lower, if the cut is relaxed into the charm region). This high-energy part of the spectrum has 
already been measured with good precision. 

Previous authors [2, 8, 9, 10, 11] have considered relations such as (1) in the slightly 
different form 




exp. input theory 

Normalizing the photon spectrum by the total 1 rate T^E*) as done in (1) has several advan- 
tages. Firstly, it is a known fact that event fractions in B — > X s ^ decay can be calculated with 
better accuracy than partial decay rates (see [13] for a recent discussion), and likewise the nor- 
malized rate does not suffer from the relatively large experimental error on the total branching 
ratio. Secondly, relation (1) is independent of the CKM factor |V#,V^*|. Thirdly, unlike the 

1 Due to an unphysical soft-photon singularity the total decay rate is commonly defined to include all events 
with photon energies above £7* = ra^j '20 [12]. 
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total B — > X s 7 decay rate, the shape of the photon spectrum is rather insensitive to possible 
New Physics contributions [12], which could distort the outcome of a \ V u b\ measurement via re- 
lation (2). Lastly, as we will see below, the weight function W(A, P + ) possesses a much better 
perturbative expansion than the function W(A,P + ) = \VtbV t *\ 2 W(A, P + )/r s (E*). This last 
point can be traced back to the fact that most of the very large contribution from the 0\ — 0-j 1 
operator mixing in the effective weak Hamiltonian cancels in the theoretical expression for the 
normalized photon spectrum. 

In principle, any partial B — > X u l~v decay rate can be brought into the form (1), with 
complicated weight functions. The relation between the two P + spectra is particularly simple, 
because the leading-power weight function is a constant at tree level. Experiments typically 
reject semileptonic 5-decay events with very low lepton energy. The effect of such an additional 
cut can be determined from [7]. Alternatively, it is possible to modify the weight function so 
as to account for a lepton cut, however at the expense of a significant increase in complexity. 
We will not pursue this option in the present work. 

The weight function depends on the kinematical variable P + and on the size A of the in- 
tegration domain. It possesses integrable singularities of the form «"(//) \n k [m b (A — P + )//i 2 ], 
with k < n, in perturbation theory. Different strategies can be found in the literature con- 
cerning these logarithms. Leibovich et al. resummed them by identifying /i 2 with m^A — P + ) 
[8, 9]. The P + dependence of the weight function then enters via the running coupling 
a s (y/mi)(A — P+j). This is a legitimate choice of scale as long as (A — P + ) has a generic value 
of order Aq CD ; however, it is not a valid choice in the small region where m b (A — P + ) ~ Aq CD . 
A key result underlying relation (1) is that, by construction, the weight function is insensitive 
to soft physics. Quark-hadron duality ensures that the region near the point P + = A, which 
is without any physical significance, does not require special consideration after integration 
over P + . In the approach of [8, 9], the attempt to resum the above logarithms near the end- 
point of the P + integral leads to integrals over unphysical Landau singularities of the running 
coupling in the nonperturbative domain. Hoang et al. chose to calculate the weight function 
in fixed-order perturbation theory at the scale /x = mj, [11]. This leads to parametrically large 
logarithms, since A— P + <C m b . In the present paper we separate physics effects from two para- 
metrically distinct scales, a hard scale Hh ~ m& and an intermediate scale /i« ~ a/^&Aqcd, so 
that we neither encounter Landau singularities nor introduce parametrically large logarithms. 
The shape of the weight function is then governed by a perturbative expansion at the inter- 
mediate scale, \n h [mi)(A — P + )//i 2 ]. As will be explained later, the coefficients in this 
series, as well as the overall normalization, possess themselves an expansion in a s (nh)- 

The calculation of the weight function starts with the theoretical expressions for the P + 
spectra in B — > X u l~v and B — > X s ^y decays, which are given as [7] 

rir r 2 it/ I 2 

x J\yy 2 - 2a ^^[(3 - 2y) ^(P + ,y) + 6(1 - y) ^ 2 (P + ,y) + y ? 3 (P+,y)} , 
= aG2 f2^ S? ^'ri {Mb ~ [Cf^ h )} 2 ^(P + ) , (4) 
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where 

P--P+ ,~ 

with P_ = Ex + | -Fx | is a kinematical variable that is integrated over the available phase 
space. Expressions for the structure functions 9^ valid at next-to-leading order (NLO) in 
renormalizat ion-group (RG) improved perturbation theory and including first- and second- 
order power corrections can be found in [7] (see also [13, 14, 15]). Symbolically, they are 
written as H(jih) • J (fa) <S> S(fa), where H(fj,h) contains matching corrections at the hard 
scale Hh- The jet function J (fa), which is a perturbative quantity at the intermediate scale 
fii, is convoluted with a non-perturbative shape function renormalized at that same scale. 
Separation of the two scales fih and /_tj allows for the logarithms in matching corrections to 
be small, while logarithms of the form ln// ft ///j, which appear at every order in perturbation 
theory, are resummed in a systematic fashion and give rise to the RG evolution functions 
U(fx h ,fa) and a r (fa, fa) [15]. 

The leading-power jet function J(p 2 , fa) entering the expressions for 3^ is universal and has 
been computed at one-loop order in [14, 15]. More recently, the two-loop expression for J has 
been obtained apart from a single unknown constant [16], which is the two-loop coefficient of 
the local S(p 2 ) term. This constant does not enter in the two-loop result for the weight function 
W(A,P + ) in (1). Due to the universality of the leading-power jet function, it is possible to 
calculate the complete 0(ot 2 s (fa)) corrections to the weight function. However, the extraction 
of hard corrections at two-loop order would require multi-loop calculations for both decay 
processes, which are unavailable at present. As a result, we will be able to predict the A and 
P + dependence of the weight function W(A, P + ) at next-to-next-to-leading order (NNLO) 
in RG-improved perturbation theory, including exact two-loop matching contributions and 
three-loop running effects. However, the overall normalization of the weight function will have 
an uncertainty of 0(a 2 (fj,h)) from yet unknown hard matching corrections. 

The total rate T S (E*) has been calculated in a local operator product expansion and reads 
(including only the leading non-perturbative corrections) [13] 



r s (£*) = 32^4 ™lm 2 M [c«(jm h )] 2 \hm\ 2 hm 



l - 



A 2 Ci(^ h ) 



9m 2 c CfJ^h) 



(6) 



where H s is the hard function of £F 7 , and H-p contains the remaining radiative corrections. 
We will present an explicit expression for this quantity at the end of Section 3 below. The 
hadronic correction proportional to A 2 /m^ cancels against an identical term in 9 r 7 . Apart 
from two powers of the running 6-quark mass defined in the MS scheme, which is part of the 
electromagnetic dipole operator O-j^ in the effective weak Hamiltonian, three more powers of 
mb emerge from phase-space integrations. To avoid the renormalon ambiguities of the pole 
scheme we use a low-scale subtracted quark-mass definition for m b . Specifically, we adopt the 
shape-function mass mf F (fa, fa) [15, 17] defined at a subtraction scale fa = 1.5 GeV, which 
relates to the pole mass as 

pole SF/ \ I CFOts(fai) ( „\ 

m? b = ml {fa, fa) +fa h.... (7) 

71 
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Throughout this paper we will use mf F ((j,*, as the 5-quark mass and refer to it as m h for 
brevity. The present value of this parameter is m b = (4.61 ± 0.06) GeV [16]. 



2 Calculation of the weight function 
2.1 Leading power 

The key strategy for the calculation of the weight function is to make use of QCD factorization 
theorems for the decay distributions on both sides of (1) and to arrange the resulting, factorized 
expressions such that they are both given as integrals over the shape function, J dui S{ui) gi(uj), 
with different functions & for the left-hand and right-hand sides. Relation (1) can then be 
enforced by matching g LH s to (?rhs- Following this procedure, we find for the integrated P + 
spectrum in B — > X u l~v decay after a series of integration interchanges 

A 1 A P+ 

JdP + ^K J dyy~ 2a H u (y) J dP + (M B - P + ) 5 J duym b J(ym b (P + - u>)) S(u) 



A M B —u> 1 

= J dwS{io) j dqhq' jdyy~ 2a H u (y)j (in m » ( ^ ? ~ ^ + Iny) , (8) 



where A g = min(A, M B — q), and [16] 

J ( ln /J'^) = / dp 2 J (p 2 , fii) (9) 

is the integral over the jet function. For the sake of transparency, we often suppress the explicit 
dependence on ^ and fih when it is clear at which scales the relevant quantities are defined. 
The function H u is a linear combination of the hard functions entering the structures 5Fj in 
(3), which in the notation of [7] is given by 

H u (y, fjL h ) = 2y 2 (3 - 2y) H ul {y, fi h ) + 12y 2 {l - y) H u2 {y, fi h ) + 2y 3 H u3 {y, fi h ) . (10) 

RG resummation effects build up the factor y~ 2a in (8), where a = ar(/J>h, Hi) is the value of 
the RG-evolution function 

arifih, to) = — r cusp (« s ( / u)) = - / da c " sp , (11) 



cusp \,"-s \t" > > I at 

which depends only on the cusp anomalous dimension [18, 19]. The quantity a has its origin 
in the geometry of time-like and light-like Wilson lines underlying the kinematics of inclusive 
B decays into light particles. Our definition is such that a is a positive number for fih > Hi 
and vanishes in the limit /i^ — ■> /ij. We find it convenient to treat the function ar((J>h, Hi) as a 
running "physical" quantity, much like a s (fi) or fn b (/j,). Since the cusp anomalous dimension 
is known to three- loop order [20] , the value of a can be determined very accurately. Note that 
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three-loop accuracy in a (as well as in the running coupling a s ) is required for a consistent 
calculation of the weight function at NNLO. The corresponding expression is 



+ 



In 



a a (fi h ) V r o A>/ 4?r 



327T 2 



r 2 


fa 




( T1 - 


A 


r 






lr 


/So 



+ . . . 



(12) 



where the expansion coefficients T n and /3 n of the cusp anomalous dimension and /5-function 
can be found, e.g., in [13]. 

Instead of the jet function J itself, we need its integral j (In Q 2 / fi 2 , //j) in the second line of 
(8). Since the jet function has a perturbative expansion in terms of "star distributions" , which 
are logarithmically sensitive to the upper limit of integration [21], it follows that j(L,/ij) is 
a simple polynomial in L at each order in perturbation theory. The two-loop result for this 
quantity has recently been computed by solving the integro-differential evolution equation 
for the jet function [16]. An unknown integration constant of 0(a 2 L°) does not enter the 
expression for the weight function. 

We now turn to the right-hand side of (1) and follow the same steps that lead to (8). 
It is helpful to make an ansatz for the leading-power contribution to the weight function, 
W^(A, P+), where the dependence on A is solely given via an upper limit of integration. To 
this end, we define a function f(k) through 



W {0 \A,P + ) oc 



1 



/•A-P+ 
JO 



dkf(k) (M B -P+-k) 1 



This allows us to express the weighted integral over the B — > X s ^ photon spectrum as 

A A P+ 

JdP + ^W^(A,P + ) oc JdP + (M B -P + ) 3 W^(A,P + ) Jdum b J(m h {P + -u;))S(Lb) 



(13) 



M 



B -u> 



oc / du S(u) 



d q W I dkf(k)j[ln mb{A « ^ kY ). (14) 



Note that the jet function J (and with it j) is the same in semileptonic and radiative decays. 
The difference is that the argument of the jet function in (8) contains an extra factor of y, 
which is absent in (14). Comparing these two relations leads us to the matching condition 



I 



o V fa 



(15) 



which holds to all orders in perturbation theory and allows for the calculation of W^(A, P + ) 
via (13). The main feature of this important relation is that the particular value of Q is irrel- 
evant for the determination of f(k). It follows that, as was the case for the jet function J, the 
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perturbative expansion of f(k) in a s (/ii) at the intermediate scale involves star distributions, 
and W^(A : P + ) depends logarithmically on (A — P + ). At two-loop order it suffices to make 
the ansatz 



f(k) oc S(k) + C F 



An 



An 



\Mi/ m b]' 



(16) 



+ ... , 



where the star distributions have the following effect when integrated with some smooth 
function (f>(k) over an interval Q: 



dk 







k 



2 ' 



L" dk G 



In 



m = [ 

Jo 



(17) 



k 



/4 



A sensitivity to the hard scale jih enters into f(k) via the appearance of H u (y,/ih) in (15). 
Because of the polynomial nature of j(L, /ij), all we ever need are moments of the hard function 
with respect to lny. We thus define the master integrals 



T n (a,fi h )= / dyy 2a H u (y,fi h ) ln n y; h n (a, fi h ) 
Jo 



T n {a,Hh) 
To(a,n h ) ' 



(18) 



which can be calculated order by order in a s (fih)- Therefore, the coefficients cj^ of the 
perturbative expansion in (16) at the intermediate scale have the (somewhat unusual) feature 
that they possess themselves an expansion in a s (/ih)- This is a consequence of the fact that, 
unlike the differential decay rates (3) and (4), the weight function itself does not obey a simple 
factorization formula, in which the hard correction can be factored out. Rather, as can be 
seen from (15), it is a convolution of the type W = H{fi h ) ® J (fa). To one-loop accuracy, the 
hard function H„ reads 



H u (y,fa) = 2y 2 (3-2y) 



1 + 



C F a s (fa 
An 



41n 2 ^ + 101n^-41n 2/ 



fa 



n 



-4L 2 (l-y)- — -12 



fa 



C F Ot s (fa) ., _> 



n 



3y lny, 



(19) 



Explicit expressions for the quantities T , c£\ and h n entering the distribution function f(k) 
will be given below. 
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2.2 Sub leading power 

Power corrections to the weight function can be extracted from the corresponding contributions 
to the two P + spectra in (3) and (4). There exists a class of power corrections associated with 
the phase-space prefactors {Mb — P+) n in these relations, whose effects are treated exactly in 
our approach, see e.g. (13). This is important, because these phase-space corrections increase 
in magnitude as the kinematical range A over which the two spectra are integrated is enlarged. 
One wants to make A as large as experimentally possible so as to increase statistics and justify 
the assumption of quark-hadron duality, which underlies the theory of inclusive B decays. 

The remaining power corrections fall into two distinct classes: kinematical corrections that 
start at order a s and come with the leading shape function [12, 21], and hadronic power 
corrections that start at tree level and involve new, subleading shape functions [23, 24, 25, 26, 
27, 28, 29, 30]. Because different combinations of these hadronic functions enter in B — > X u l~v 
and B — > X s ^ decays, it is impossible to eliminate their contributions in relations such as 
(1). As a result, at 0(AQcr>/ m b) there are non-perturbative hadronic uncertainties in the 
calculation of the weight function W(A, P+), which need to be estimated before a reliable 
extraction of \V u t>\ can be performed. For the case of the charged-lepton energy spectrum and 
the hadronic invariant mass spectrum, this aspect has been discussed previously in [25, 26] 
and [27], respectively. 

Below, we will include power corrections to first order in Aqcd /™>b- Schematically, the 
subleading corrections to the right-hand side of (1) are computed according to ~ ® 
dT^P /dP + + <S>dT^ /dP + , where the superscripts indicate the order in 1/mj, power count- 
ing. The power corrections to the weight function, denoted by are derived from the 
mismatch in the power corrections to the two decay spectra. The kinematical power correc- 
tions to the two spectra are known at 0(a s ), without scale separation. We assign a coupling 
a s (fl) to these terms, where the scale Jx will be chosen of order the intermediate scale [7]. At 
first subleading power the leading shape function is convoluted with either a constant or a 
single logarithm of the form ln[(P + — Cj)/{Mb — P+)], and we have (with n — 0, 1) 



/ dP + —^3 a s (fi) / dP + {M B - P+) 4 / du S(u) In" 
Jo dr+ Jo Jo 



M B -P 4 



= a s (p) dujS(uj) dk (M B - Cj - kf In™ -, (20) 

Jo Jo M B -uj-k 

and similarly for the photon spectrum. On the other hand, the weighted integral in (1) also 
contains terms where the photon spectrum is of leading power and the weight function of 
subleading power, 

/ dP + --^W kin( - 1) {A,P + ) 3 / dP + {M B -P + ) 3 S{P + )W kin(1 \A,P + ). (21) 
Jo dP+ Jo 

Therefore the kinematical corrections to the weight function must have the form 

w ^ a - p +> °= sr ,,k (mb - p+ - kr ( A+Bia M B - k p+-k ) • <22> 
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and a straightforward calculation determines the coefficients A and B. The hadronic power 
corrections to the weig ht function, W hadl ^\ can be expressed in terms of the subleading shape 
functions t(u), u(ou), and v(u) defined in [29]. These terms are known at tree level only, and 
at this order their contribution to the weight function can be derived using the results of [7]. 



3 Results 



Including the first-order power corrections and the exact phase-space factors, the weight func- 
tion takes the form 



W(A,P + ) 



G 2 F ml 

1927T 3 



T {a,n h )H T (n h ) (M B -P + 



x 



1 + 



47T 







+ c 



In 



m b (A - P + ) 



47T 



+ 4 2) ( In 



m b (A - P + ) 



2 Pl {5) In 



m 6 (A - P + ) 



+ 2p 2 (S) 



+ , A ^ 9l^± [ A (a, » h ) I A (S) + B(a, Ib{5)\ 



Mr - P. 4tt 



+ 



M B - P + 2(l-a)(3-a) 



4(1 -a)(A-P+) + 2(4-3a) 



i(P + ) 
S(P + ) 



+ (4 



S(P + ) 



+ (8 - 13a + 4a 2 



HP + ) 
S(P + ) 



+ 



Mn-P+ 



s\p+) 

S{P + ) 



(23) 



+ . . . 



where A = Mb — ray, is the familiar mass parameter of heavy-quark effective theory, and 
5 = (A — P + )/ (M B — P+). The first line denotes an overall normalization, the next three lines 
contain the leading-power contributions, and the remaining expressions enter at subleading 
power. The different terms in this result will be discussed in the remainder of this section. 

For the leading-power terms in the above result we have accomplished a complete sep- 
aration of hard and intermediate (hard-collinear) contributions to the weight function in a 
way consistent with the factorization formula W = H <8> J mentioned in the previous section. 
The universality of the shape function, which encodes the soft physics in both B — > X s ^ and 
B — > X u l~v decays, implies that the weight function is insensitive to physics below the in- 
termediate scale ^ ~ y/m&AQCD- In particular, quark-hadron duality ensures that the small 
region in phase space where the argument ra^/S. — P + ) of the logarithms scales as Aq CD or 
smaller does not need special consideration. At a technical level, this can be seen by noting 
that the jet function is the discontinuity of the collinear quark propagator in soft-collinear 
effective theory [15, 22], and so the P + integrals can be rewritten as a contour integral in the 
complex p 2 plane along a circle of radius m& A ~ fif. 
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3.1 Leading power 



The leading-power corrections in the curly brackets in (23) are determined completely at 
NNLO in RG-improved perturbation theory, including three-loop running effects via the quan- 
tity a in (fl), and two-loop matching corrections at the scale fii as indicated above. To this 
end we need expressions for the one-loop coefficients c& including terms of 0(a s (j2f l )), while 

(2) 

the two-loop coefficients c„ are needed at leading order only. We find 



-3hi(a,fx h ) 



1 - 



C F a s ([i h ) 4/i* 



7T 



3m& 



+ 2h 2 (a,Hh) , 



^(a,^), (24) 



and 



J2) 



= c, 



- + 27r 2 -24C 3 hi^fi^+i - 



9 4tT 



.(2) 



+ C A 

+ Po 
C F 



73 

-— +40C 3 ) Ma,W») + 



2 3 
8 2tt 2 



h 2 (a, n h ) - 6h 3 (a, fx h ) + 2h 4 (a, fi h ) 



h 2 {a,fj, h ) 



247 2tt 2 
+ 



18 3 
12h 2 (a,n h ) + Sh 3 (a,fj, h ) 
29 

— hi(a,fj, h ) - 2h 2 (a,n h ) 



29 2 
hi(a, n h ) + — h 2 {a, fi h ) - - h 3 {a, fji h ) 



3 

16 4tt 2 



hi(a,/j, h ) 



.(2) 



8C F h 2 (a, fx h ) - 2(3o h^a, fj, h ) . 



(25) 



As always C F = 4/3, Ca = 3, and [3 = 11— 2n//3 is the first coefficient of the QCD /9-function. 
The term proportional to /i* in the expression for arises because of the elimination of the 
pole mass in favor of the shape-function mass, see (7). Since the logarithms ln[m&(A — P + )/fif] 
in (23) contain m&, all coefficients except receive such contributions. However, to two-loop 
order only is affected. 

Next, the corresponding expressions for the hard matching coefficients hi are calculated 
from (18). To the required order they read 



hi(a,fj, h ) 



15- 12a + 2a 2 



2(2-a)(3-a)(3-2a) 



+ 



C F a s (fi h ) 
An 



2(189 - 318a + 192a 2 - 48a 3 + 4a 4 ) , m b 

— in — 



(2-a) 2 (3-a) 2 (3-2a)^ 



2331 - 5844a + 5849a 2 - 2919a 3 + 726a 4 - 72a 5 ,, w 
+ (2- B ).(3- B H3-2 B ). 4^)(3-2a) 



+ 



9 



h 2 {a,fj, h ) 



69 - 90a + 36a 2 - 4a 3 
2(2-a) 2 (3-a)(3-2a) 2 

C F a s (n h ) 



+ 



Att 



2(1692 - 3699a + 3138a 2 - 1272a 3 + 240a 4 - 16a 5 ) m b 
(2-a) 3 (3-a) 2 (3-2a) 3 H 7^ 



46521 - 140064a + 175479a 2 - 117026a 3 + 43788a 4 - 8712a 5 + 720a 6 

2(2-a) 4 (3-a) 2 (3-2a) 4 

4(15 - 12a + 2a 2 ) , W/ . , 
+ (2- a )(3- a )(3-L) * 1(3 " 2 "» " 4 *' >(3 " 2a) + ' ' ' 

3(303 - 552a + 360a 2 - 96a 3 + 8a 4 ) 
h ^^)~ 4(2-a) 3 (3-a)(3-2a) 3 + '-" 

, , , 3(1293 - 3030a + 2760a 2 - 1200a 3 + 240a 4 -16a 5 ) , , 

k ^ ^ = 2(2-a) 4 (3-a)(3-2a) 4 + " " " ' (26) 

where (x) is the n-th derivative of the polygamma function. Because of the exact treatment 
of the phase space there are corrections to the logarithms in (23), which are finite-order 
polynomials in the small ratio 5 = (A — P + )/(M B — P+). Explicitly, 

»W=W-§* + £*-^ + i*. (27) 
This concludes the discussion of the leading-power expression for the weight function. 

3.2 Sub leading power 

The procedure for obtaining the kinematical power corrections to the weight function has been 
discussed in Section 2.2. For the coefficients A and B in (23) we find 

-388 + 702a - 429a 2 + 123a 3 - 34a 4 + 8a 5 (\ 4 , m b \ [Cf?(^ ft )] : 
A ( a > Vh) = TTTi W77^ + o _ n ln 



2(l-a) 2 (2-a)(3-a)(3-2a) \3 9 mj [C^(fi h )] 2 

10 ggW 8 / CM 1 CMCgM \ 16 [Ci(/x h )] 2 

~ 3 G*M + 3 3 [CfiW]* ) 9l[Z) 9 [C|^)] 2 92[Z) ' 

m „„x 2(8 + a) 2 [CfMf 

(a ' ^ " "(l-a)(3-a) " 9 \C^W ' ( §) 

Here Cj(/i/i) denote the (effective) Wilson coefficients of the relevant operators in the effective 
weak Hamiltonian, which are real functions in the Standard Model. The variable z = (m c /mi>) 2 
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enters via charm-loop penguin contributions to the hard function of the B — > X s ^ photon 
spectrum [12], and 



9iW 



dxxRe 



G 



x 



(!) + 5 



dx (1 — x ) 



G 



x 



+ 5 



with 



-2arctan 2 - N /t/(4 -t) 



G(t) = 



2 In 



■ t <4, 

• \ 2 

y 



(29) 



(30) 



Furthermore we need the integrals over k in (22), which encode the phase-space corrections. 
They give rise to the functions 



I A {8) = l-25 + 25 2 -o 3 + -o 4 , 

5 

T /n t /r\ i 5 ln(l -5) 4 3 . 4 r2 1 3 

J B (<$) = J A (<$) In H ^— — \--S (TH 5 3 . 

BK ' Ay ' 1-5 55 5 5 15 20 



(31) 



The hadronic power corrections come from subleading shape functions in the theoretical 
expressions for the two decay rates. We give their tree-level contributions to the weight 
function in the last two lines of (23), where S denotes the leading shape function, and t, u, v are 
subleading shape functions as defined in [29]. For completeness, we also include a contribution 
proportional to m 2 s resulting from finite-mass effects in the strange-quark propagator in B — > 
X s 7 decays. For m s = 0(Aq C d) these effects are formally of the same order as other subleading 
shape-function contributions [31], although numerically they are strongly suppressed. The 
appearance of subleading shape functions introduces an irreducible hadronic uncertainty to a 
\V u b\ determination via (1). In practice, this uncertainty can be estimated by adopting different 
models for the subleading shape functions. This will be discussed in detail in Section 4.4 
below. Until then, let us use a "default model" , in which we assume the functional forms of 
the subleading shape functions t(ui), w(cD), and v(u) to be particular linear combinations of 
the functions S'(u) and (A — Co) S(Cu). These combinations are chosen in such a way that the 
results satisfy the moment relations derived in [29], and that all terms involving the parameter 
A cancel in the expression (23) for the weight function for any value of a. These requirements 
yield 



t 



- (A 

4 V 



A 2 + — 



2 v ; 6 



A 2 5", (32) 



and the last two lines inside the large bracket in the expression (23) for the weight function 
simplify to 

M) S'(P + ) , Al SF (a) un a \ 4A! SF (a) 



where 



M B ~ P+ S{P+) M B - P_ 
(2 + a) Ai 



5(P + - A) 



(M B - P+) 



2 ' 



a ssf(«) 



a(7 — 4a) \ 2 
3(l-a)(3-a) ' 2(1 -a)(3-a) 



+ 



(33) 
(34) 
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Here Ai and A 2 = \{Mg» — M|) are hadronic parameters describing certain B- meson matrix 
elements in heavy-quark effective theory [32]. The strange-quark mass is a running mass 
evaluated at a scale typical for the final-state hadronic jet, for which we take 1.5 GeV. As 
mentioned above, the numerical effect of the strange-quark mass correction is small. For 
typical values of the parameters, it reduces the result for A| SF by about 10% or less. The 
expression on the right-hand side in (33) is equivalent to that on the left-hand side after the 
integration with the photon spectrum in (1) has been performed. It has been derived using 
the fact that the normalized photon spectrum is proportional to the shape function S(P + ) at 
leading order. Note that the second term in the final formula is power suppressed with respect 
to the first one. It results from our exact treatment of phase-space factors and thus is kept 
for consistency. 



3.3 Normalization 

Finally, let us present explicit formulae for the overall normalization factor in (23). The new 
ingredient here is the factor T , which is defined in (18). At one-loop order we find 



T (a,n h ) 



2(3 - a) 



(2-a)(3-2a) 



1 - 



Att 



4 In 



2 m b 2(120 - 159a + 69a 2 - 10a 3 ) m b 
(2 - a)(3 - a)(3 - 2a) 



7T 



1539 - 3845a + 3842a 2 - 1920a 3 + 480a 4 - 48a 5 , , m . 
+ (2- a )»(3-a)(3-2.)» + ^ >(3 " 2a) + 6 



When the product of T with the quantity [13] 

C F a s (fx h 



HrM = 1 + 



47T 



41n 2 ^-101n^ + 7 



7tt 2 | 12/i* 



6 rrib 
2 



2 In 2 8, - (7 + AS, - 5 2 ) In 6* + 105* + 5 2 - - <5 3 
+ t^-5t, — — /n(<>*) + T^ffT — r hi\P*) H FT^fTT — xv> — h*\P*) 



+ 



i fX \ , l^tki i^h)} 2 t /j- s 



(35) 



(36) 



from the total B — > X s 7 decay rate is consistently expanded to 0(a s (fih)), the double loga- 
rithm cancels out. Here 5* = 1 — 2£'*/m& = 0.9, and the functions fij(6*) capture effects from 
operator mixing. 
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Figure 1: Residual scale dependence of T M (A) for A = 0.65 GeV. The dashed line depicts the 
variation of fih about its default value m^/v^, the dash-dotted line the variation of /i« about 
1.5 GeV, and the dotted line the variation of p, also about 1.5 GeV. The highlighted area shows 
the combined perturbative uncertainty. 

4 Numerical results 

We are now in a position to explore the phenomenological implications of our results. We 
need as inputs the heavy-quark parameters A 2 = 0.12 GeV 2 , Ai = (—0.25 ± 0.10) GeV 2 , and 
the quark masses m b = (4.61 ± 0.06) GeV [16], m s = (90 ± 25) MeV [33, 34], and m c /m h = 
0.222±0.027 [13]. Here m b is defined in the shape-function scheme at a scale /i* = 1.5 GeV, m s 
is the running mass in the MS scheme evaluated at 1.5 GeV, and m c jm b is a scale invariant 
ratio of running masses. Throughout, we use the 3-loop running coupling normalized to 
a s (Mz) = 0.1187, matched to a 4-havor theory at 4.25 GeV. For the matching scales, we pick 
the default values /i^ ef = m b /^/2 and [if ef = /2 def = 1.5 GeV, which are motivated by the 
underlying dynamics of inclusive processes in the shape-function region [7, 15]. 

In the remainder of this section we present results for the partial decay rate T U (A) com- 
puted by evaluating the right-hand side of relation (1). This is more informative than to 
focus on the value of the weight function for a particular choice of P + . For the purpose of 
our discussion we use a simple model for the normalized photon spectrum that describes the 
experimental data reasonably well, namely 

with A = 0.77 GeV and b = 2.5. 

4.1 Studies of the perturbative expansion 

The purpose of this section is to investigate the individual contributions to T U (A) that result 
from the corresponding terms in the weight function, as well as their residual dependence on 
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Figure 2: Convergence of the perturbative series and residual scale dependence. Left: RG- 
improved results at LO (dotted), NLO (dashed), and NNLO (solid) as a function of the scales 
Hi, Hh, H, which are varied simultaneously about their default values. Right: Fixed-order 
results at tree-level (dotted) and one- loop order (dashed). 



the matching scales. For A = 0.65 GeV we find numerically 
r u (0.65GeV) 

IKbPps- 1 



43.5 (1 + 0.158 M«)] - 0.095 M/^)] + 0.076 

- 0.037 M^KO^)] + 0.009 [kin] - 0.043 [hadr]) = 46.5 . 



(38) 



The terms in parenthesis correspond to the contributions to the weight function arising at 
different orders in perturbation theory and in the l/m^ expansion, as indicated by the sub- 
scripts. Note that the perturbative contributions from the intermediate scale are typically 
twice as large as the ones from the hard scale, which is also the naive expectation. Indeed, the 
two- loop a 2 s (ni) correction is numerically of comparable size to the one-loop a s (Hh) contribu- 
tion. This confirms the importance of separating the scales Hi and Hh- The contributions from 
kinematical and hadronic power corrections turn out to be numerically small, comparable to 
the two-loop corrections. 

The weight function (23) is formally independent of the matching scales Hh, Hi, and p,. 
In Figure 1 we plot the residual scale dependence resulting from the truncation of the per- 
turbative series. Each of the three scales is varied independently by a factor between l/y/2 
and y/2 about its default value. The scale variation of Hi i s still as significant as the varia- 
tion of Hh, even though the former is known at NNLO and the latter only at NLO. We have 
checked analytically that the result (23) is independent of Hi through two-loop order, i.e. the 
residual scale dependence is an 0{aP s {Hi)) effect. In order to obtain a conservative estimate 
of the perturbative uncertainty in our predictions we add the individual scale dependencies in 
quadrature. This gives the gray band shown in the figure. 

Figure 2 displays the result for r„(0.65GeV) at different orders in RG-improved pertur- 
bation theory. At LO, we dismiss all a s terms including the kinematical power corrections; 
however, leading logarithms are still resummed and give rise to a non-trivial dependence of 
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T on the coefficient a. At NLO, we include the 0(a s (fj,f l )), 0(a s (/ij)), and 0(a s (p,)) con- 
tributions, but drop terms of order ot 2 s {[ii) or ct s (//j) a s {yi,h}. At NNLO, we include all terms 
shown in (23). In studying the different perturbative approximations we vary the matching 
scales simultaneously (and in a correlated way) about their default values. Compared with 
Figure 1 this leads to a reduced scale variation. The gray bands in Figure 2 show the total 
perturbative uncertainty as determined above. While the two-loop NNLO contributions are 
sizable, we observe a good convergence of the perturbative expansion and a reduction of the 
scale sensitivity in higher orders. The right-hand plot in the figure contrasts these findings 
with the corresponding results in fixed-order perturbation theory, which are obtained from 
(23) by setting fj, h — ^ — p, — /i and truncating the series at 0(a s (fi)) for consistency. We 
see that the fixed-order results are also rather insensitive to the value of /i unless this scale 
is chosen to be small; yet, the predicted values for T u are significantly below those obtained 
in RG-improved perturbation theory. We conclude that the small scale dependence observed 
in the fixed-order calculation does not provide a reliable estimator of the true perturbative 
uncertainty. In our opinion, a fixed-order calculation at a high scale is not only inappropriate 
in terms of the underlying dynamics of inclusive decay processes in the shape-function region, 
it is also misleading as a basis for estimating higher-order terms in the perturbative expansion. 

4.2 Comments on the normalization of the photon spectrum 

We mentioned in the Introduction that the use of the normalized photon spectrum is ad- 
vantageous because event fractions in B — > X s 7 decay can be calculated more reliably than 
partial decay rates. In this section we point out another important advantage, namely that 
the perturbative series for the weight function W(A, P + ) is much better behaved than that 
for W(A,P + ). The difference of the two weight functions lies in their normalizations, which 
are 



Here H s is the hard function in the factorized expression for the structure function 5F 7 in 
(4), which has been derived in [13]. Note that the two weight functions have a different 
dependence on the 6-quark mass. In the case of W, three powers of enter through phase- 
space integrations in the total decay rate T S (E*), and it is therefore appropriate to use a low- 
scale subtracted quark-mass definition, such as the shape-function mass. In the case of W, on 
the other hand, two powers of the running quark mass m^/x/J enter through the definition of 
the dipole operator 0-j 1 ^ and it is appropriate to use a short-distance mass definition such as 
that provided by the MS scheme. In practice, we write mb(Hh) as rfibirrib) times a perturbative 
series in a s (/ih)- 

The most pronounced effect of the difference in normalization is that the weight function 
W receives very large radiative corrections at order a s (nh), which range between —68% and 
—43% when the scale Hh is varied between m;, and rrib/2. This contrasts the well-behaved 
perturbative expansion of the weight function W, for which the corresponding corrections 
vary between —11% and —7%. In other words, the hard matching corrections for W are about 



W{A, P+) oc ml T (a, // h )#rM , W(A, P + ) oc 



T (a,/j, h ) 



(39) 
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Figure 3: Perturbative uncertainties on T n (A) encountered when using the weight function 
W{ A, P + ) (dark gray) or W(A,P + ) (light gray). In the latter case, the normalization of the 
photon spectrum is chosen such that the two predictions agree at A = 0.65 GeV for central 
values of the matching scales. 



six times larger than those for W. Indeed, these corrections are so large that in our opinion 
relation (2) should not be used for phenomenological purposes. 

The different perturbative behavior of the hard matching corrections to the weight func- 
tions is mostly due to the mixing of the dipole operator Oy 7 with other operators in the 
effective weak Hamiltonian for B — > X s ^ decay. In order to illustrate this fact, consider the 
one-loop hard matching coefficients defined as 

^ocl + fc^U..., Wocl + k^ + .... (40) 

7T 7T 

With our default scale choices we have k = —2.32 + 1.13 = —1.19, where the second contribu- 
tion (+1.13) comes from operator mixing, which gives rise to the terms in the last two lines 
in (36). For the weight function W, this contribution has the opposite sign than the other 
terms, so that the combined value of k is rather small. For the weight function W, on the 
other hand, we find k = —1.58 — 5.52 = —7.10. Here the contribution from operator mixing is 
dominant and has the same sign as the remaining terms, thus yielding a very large value of k. 
Such a large 0(a s ) correction was not observed in [11], because these authors chose to omit 
the contribution from operators mixing. Note that at a higher scale Hh = Tn h , as was adopted 
in this reference, the situation is even worse. In that case we find k = —2.81 + 1.19 = —1.62 
and k = -0.66 - 9.13 = -9.79. 

A visualization of the perturbative uncertainty is depicted in Figure 3, where predictions 
for r u (A) are shown using either (1) or (2). In each case, the error band is obtained by 
varying the different scales about their default values, /i„ G [fj,^ /\/2, /x£ ef \/2], and adding the 
resulting uncertainties in quadrature. The dark-gray band bordered by solid lines denotes 
the perturbative uncertainty of predictions when using the normalized photon spectrum, as 
in (1). (At the point A = 0.65 GeV this uncertainty is identical to the gray band depicted 
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in Figures 1 and 2). The light-gray band bordered by dashed lines corresponds to the use of 
the absolute photon spectrum, as in (2). The difference in precision between the two methods 
would be even more pronounced if we used the higher default value fih — n^b f° r the hard 
matching scale. Obviously, the use of the normalized photon spectrum will result in a more 
precise determination of |14&|- 



4.3 Comments on (3qOl 2 terms and scale separation 



s 



The separation of different momentum scales using RG techniques, which is one of the key 
ingredients of our approach, is well motivated by the dynamics of charmless inclusive B decays 
in the shape-function region. Factorized expressions for the -B-decay spectra involve hard 
functions renormalized at fih multiplied by jet and shape functions defined at a lower scale. 
While physics at or below the intermediate scale is very similar for B — > X u l~v and B — > 
X s 7 (as is manifested by the fact that the leading shape and jet functions are universal), 
the physics at the hard scale in B — > X s ^ decay is considerably more complicated than in 
semileptonic decay, and it might even contain effects of New Physics. Therefore it is natural 
to respect the hierarchy Hh ^> //j and disentangle the various contributions, as done in the 
present work. In fact, our ability to calculate the dominant two- loop corrections is a direct 
result of this scale separation. Nevertheless, at a technical level we can reproduce the results 
of a fixed-order calculation by simply setting all matching scales equal to a common scale, 

= (J>i = p> = f-i- In this limit, the expressions derived in this work smoothly reduce to those 
obtained in conventional perturbation theory. While factorized expressions for the decay rates 
are superior to fixed-order results whenever there are widely separated scales in the problem, 
they remain valid in the limit where the different scales become of the same order. _ 

In a recent publication, the O(f3oa 2 ) BLM corrections [35] to the weight function W(A, P + ) 
in (2) were calculated in fixed-order perturbation theory [11]. For simplicity, only the contri- 
bution of the operator to the B — > X s 7 decay rate was included in this work. We note that 
without the contributions from other operators the expression for W is not renormalization- 
scale and -scheme invariant. Neglecting operator mixing in the calculation of W is therefore 
not a theoretically consistent approximation. However, having calculated the exact NNLO 
corrections at the intermediate scale allows us to examine some of the terms proportional to 
f3oa 2 s {iXj) and compare them to the findings of [11]. In this way we confirm their results for 
the coefficients multiplying the logarithms ln n [m;,(A — P + )/fi 2 ] with n — 1,2 in (23). While 
the p a 2 terms approximate the full two-loop coefficients of these logarithms arguably well, 
we stress that the two-loop constant at the intermediate scale is not dominated by terms 
proportional to flo- Numerically we find 



c {2) = -47 .4 + 39.6 + a 
25/3 



-31.8 + 38.8 '" 



25/3 



+ 0(a 2 ), (41) 



which means that the approximation of keeping only the BLM terms would overestimate 
this coefficient by almost an order of magnitude and give the wrong sign. This shows the 
importance of a complete two-loop calculation, as performed in the present work. 

We believe that the perturbative approximations adopted in our paper, i.e. working to 
NNLO at the intermediate scale and to NLO at the hard scale, are sufficient for practical 
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purposes in the sense that the residual perturbative uncertainty is smaller than other uncer- 
tainties encountered in the application of relation (1). Still, one may ask what calculations 
would be required to determine the missing a^(fih) terms in the normalization of the weight 
function in (23), or at least the terms of order (3 a^(fih)- For the case of B — > X s ^ decay, the 
contribution of the operator to the normalized photon spectrum was recently calculated 
at two-loop order [36], while the contributions from other operators are known to O(/3 a: 2 ) 
[37]. What is still needed are the two-loop corrections to the double differential (in P + and y) 
B — > X u l~D decay rate in (3). 

4.4 Subleading corrections from hadronic structures 

Due to the fact that different linear combinations of the subleading shape functions t(cu), 
u(oo), and v(u) enter the theoretical description of radiative and semileptonic decays starting 
at order Aqcd/^6, the weight function cannot be free of such hadronic structure functions. 
Consequently, we found in (23) all of the above subleading shape functions, divided by the 
leading shape function S(cu). Our default model (32) for the subleading shape functions was 
chosen such that the combined effect of all hadronic power corrections could be absorbed into 
a single hadronic parameter Ag SF . More generally, we define a function Ohadr(A) via (a factor 
2 is inserted for later convenience) 

T U (A) = [r M (A)] dcf [1 + 2o hadr (A)] , (42) 

where [r u (A)]d e f denotes the result obtained with the default model for the subleading shape 
functions. From (23), one finds that 

J Q A dP + (M B - P + Y [2(4 - 3a) h t (P + ) + (4 - a) h u (P + ) + (8 - 13a + 4a 2 ) h v (P + )\ 

Ohadr ( A ) — a , 

4(l-a)(3-a)f A dP + (M B -P + fS(P + ) 

(43) 

where we have used that, at leading order in a s and Aqcd/^6, the B — > X s ^ photon spectrum 
is proportional to (M B — P + ) 3 S(P + ). In the relation above, h t (Co) = i(uo) — [£(u))] def etc. denote 
the differences between the true subleading shape functions and the functions adopted in our 
default model. By construction, these are functions with vanishing normalization and first 
moment. 

The above expression for Ohadr(A) is exact to the order we are working; however, in practice 
we do not know the precise form of the functions hi(u). Our goal is then to find a conservative 
bound, |<5hadr(A)| < e hadr (A), and to interpret the function ehadr(A) as the relative hadronic 
uncertainty on the value of \V u b\ extracted using relation (1). To obtain the bound we scan 
over a large set of realistic models for the subleading shape functions. In [7], four different 
functions hi(u>) were suggested, which can be added or subtracted (in different combinations) 
to each of the subleading shape functions. Together, this provides a large set of different 
models for these functions. To be conservative, we pick from this set the model which leads 
to the largest value of |5hadr(A)|. The integrand in the numerator in (43) is maximized if all 
three hi(Co) functions are equal to a single function, whose choice depends on the value of A. 
In the denominator, we find it convenient to eliminate the shape function S(P + ) in favor of 
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Figure 4: Estimates for the hadronic uncertainty ehadr (A) obtained from a scan over models 
for the subleading shape functions. The dashed lines correspond to the individual results for 
the four hi(u) functions suggested in [7]. The thick solid line, which covers one of the dashed 
lines, shows the maximum effect. 



the normalized photon spectrum. Working consistently to leading order, we then obtain 



£hadr(A) — 



5 — 5a + a 2 U(/ih,fa) 



[l-a)(S-a) 



max 



dP + {M B -P + ) A h t {P+ 



mi 



f 

Jo 



dP + (M B - P + 



1 



aT s 



(44) 



r,(Ej dP+ 

where as before a = ar(fai, fa) ~ 0.12 for the default choice of matching scales, and U (fa, fa) ~ 
1.11 [7]. 

The result for the function e iiadr (A) obtained this way is shown in Figure 4. We set the 
matching scales to their default values and use the model (37) for the photon spectrum, which 
is a good enough approximation for our purposes. From this estimate it is apparent that the 
effects of subleading shape functions are negligible for large values A 3> Aqcd and moderate 
for A ~ Aqcd, which is the region of interest for the determination of \V u b\. In the region 
A < 0.3 GeV the accuracy of the calculation deteriorates. For example, we find e ha dr = 2.0% 
for A = 0.65 GeV, e hadr = 4.8% for A = 0.5 GeV, and e hadr = 8.2% for A = 0.4 GeV. 



5 Conclusions 

Model-independent relations between weighted integrals of B — > X s 7 and B — > X u l~v decay 
distributions, in which all reference to the leading non-perturbative shape function is avoided, 
offer one of the most promising avenues to a high-precision determination of the CKM matrix 
element \V u b\. In order to achieve a theoretical precision of better than 10%, it is necessary to 
include higher-order corrections in a s and Aq CD /m fe in this approach. 

In the present work, we have calculated the weight function W(A, P + ) in the relation 
between the hadronic P + = Ex — \Px\ spectra in the two processes, integrated over the 
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interval < P + < A. Based on QCD factorization theorems for the differential decay rates, 
we have derived an exact formula (15) that allows for the calculation of the leading-power 
weight function to any order in perturbation theory. We have calculated the A- and P + - 
dependent terms in the weight function exactly at next-to-next-to-leading order (NNLO) in 
renormalizat ion-group improved perturbation theory, including two-loop matching corrections 
at the intermediate scale /Zj ~ iy/m b AQ CD and three-loop running between the intermediate 
scale and the hard scale fih ~ m^. The only piece missing for a complete prediction at 
NNLO is the two-loop hard matching correction to the overall normalization of the weight 
function. A calculation of the a"i(fih) term would require the knowledge of both decay spectra 
at two-loop order, which is currently still lacking. We also include various sources of power 
corrections. Power corrections from phase-space factors are treated exactly. The remaining 
hadronic and kinematical power corrections are given to first order in Aqcd/^6 and to the 
order in perturbation theory to which they are known. 

A dedicated study of the perturbative behavior of our result for the weight function has 
been performed for the partial B — > X u l~u decay rate T M (A) as obtained from the right- 
hand side of relation (1). It exhibits good convergence of the expansion and reduced scale 
sensitivity in higher orders. We find that corrections of order a^(fii) at the intermediate scale 
are typically as important as first-order a s (fih) corrections at the hard scale. We have also 
seen that fixed-order perturbation theory significantly underestimates the value of r u (A), 
even though the apparent stability with respect to scale variations would suggest a good 
perturbative convergence. In order to obtain a well-behaved expansion in powers of a s (fih), it 
is important to use the normalized photon spectrum in relation (1). A similar relation involving 
the differential B — > X s 7 decay rate receives uncontrollably large matching corrections at the 
hard scale and is thus not suitable for phenomenological applications. At next-to-leading 
order in the 1 /nib expansion, the weight function receives terms involving non-perturbative 
subleading shape functions, which cannot be eliminated. Our current ignorance about the 
functional forms of these functions leads to a hadronic uncertainty, which we have estimated by 
scanning over a large set of models. We believe that a reasonable estimate of the corresponding 
relative uncertainty ehadr on \V u b\ is given by the solid line in Figure 4. 

Let us summarize our main result for the partial B — > X u l~u decay rate with a cut 
P+ < 0.65 GeV, which is close to the charm threshold M^/Mb, and present a detailed list of 
the various sources of theoretical uncertainties. We find 

r u (0.65 GeV) = (46.5 ± 1.4 [pert] ± 1.8 [hadr] ± 1.8 [m b ] ± 0.8 [pars] ± 2.8 [norm] ) | V^fo | 2 pS^ 1 

= (46.5 ±4.1) \V ub \ 2 ps-\ (45) 

where the central value is derived assuming that the B — > X s ^ photon spectrum can be 
accurately described by the function (37). The errors refer to the perturbative uncertainty as 
estimated in Section 4.1, the uncertainty due to the ignorance about subleading shape functions 
as discussed in Section 4.4, the error in the value of the 6-quark mass, other parametric 
uncertainties from variations of m c , m s , and Ai, and finally a 6% uncertainty in the calculation 
of the normalization of the photon spectrum [13] . To a good approximation the errors scale 
with the central value. The above numbers translate into a combined theoretical uncertainty 
of 4.4% on \V u b\ when added in quadrature. 
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